Chandra Contaminant Migration Model 

Goals: 

(1) Determine if the observed change in Optical Blocking Filet (OBF) 
accumulation rate and spatial distribution can be explained by the 
known change in thermal environment caused by the ACIS housing 
(AH) heater being turned off in 2008 April. 

(2) Constrain the material properties of the contaminant(s). 

Methods: 

(1) Simulate migration (vaporization/deposition) of potential organic 
contaminant(s) through evaporation (or sublimation) within thermal & 
geometric model constraints. 

(2) Compare time evolution of material mass column density to 
observations at select locations on OBF filters to identify best 
candidate contaminant(s). 

Study is extension of bake-out studies (2004, see O’Dell et al. SPIE 
5898,313 (2005)). Details of new work, O’Dell et al. SPIE 8859,0 
(2013); especially, improved geometric/thermal model by Neil Tice. 
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Chandra Contaminant Migration Model 

Observational Constraints: 

(1) OBF center accumulation history as determined by Alexey Vikhlinin 
(AV) from ECS and A1795 observations, expressed in terms of optical 
depth x at Mn-L (-640 eY); [AV memo Dec. 2012] 

(2) OBF edge-to-center evolution, edge at -row 64 center at row 512 
[op.cit.] 

(3) OBF near-edge (row 167) accumulation history as determined by 
Herman Marshall (HM) from grating observations, expressed in terms of 
optical depths at C and O absorption edges. Plus a few measurements 
nearer and farther from row 167. 

Contamination Migration Model reports contaminant mass column density. 
Optical depth unity at Mn-L is roughly equivalent to a 100 pg/cm 2 column 
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Chandra Contaminant Migration Model 

The basic algorithm: 

♦ Temperature dependence of vaporization rate follows Clausius- 
Clapeyron relation 

♦ p(T) = p(To) exp (-AH(To-T)/RTTo) 

where p(T) is the mass vaporization rate at temperature T, AH 

is the vaporization enthalpy, and R the universal gas constant 

♦ Solves for the evolution/migration of the mass column of 
contaminant on all surfaces using an explicit (forward) difference 
scheme to solve the initial value problem: 

♦ dfjn/dt =-p(Ti)+Zj p(Tj)(Aj/Ai)Sji 

where the Aj are areas of surfaces and the view factors Sji are the 
fraction of the opening angle of surface j subtended by surface i. 
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Chandra Contaminant Migration Model 

All simulations assume: 

♦ A source rate of the form S(t)=A+Bexp(-t/3.7) for all times, t, 
where time-dependence adjusted to follow AV’s center x evolution 
and normalization to match HM’s corresponding mass surface 
density estimates for t<9 yr. 

♦ Time-averaged surface temperatures with AH heaters ON until 
-2008.5 (t=9 yr) then AH heaters OFF as determined by Neil Tice 
(NT) from flight data & his geometric/thermal model 

n 

♦ An organic material w/ mass vaporization rate -10' 7 //g/cm 2 /s at 
temperatures —50 °C 

♦ Emissivity 8 = 0.05 to 0.40; determines OBF temperatures 
(other surfaces are conductively coupled and insensitive to e) 

♦ A fixed fraction, /; = Avent/(Avent+ATc), of contaminant striking 
the telescope closeout (TC) at the top of the geometric model 
vents to space, where A = area 
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Summary of vaporization or sublimation 


Name 

tetra- 

penta- 

hexa- 

hepta- 

octa- 

nona- 

ei- 

heni- 

do- 

tri- 

tetra- 

penta- 

tetra- 



decane 3 

decane 3 

decane b 

decane b 

decane 

decane c 

cosane 

cosane 

cosane 

cosane c 

cosane 

cosane 

cosane 

DOP 

Unit 

Formula 

C 14 H 30 

C 15 H 32 

C 16 H 34 

C 17 H 36 

00 

CO 

X 

00 

d 

C-19H40 

C20H42 

C 2 i H 44 

C22H46 

C 23 H 48 

C 24 H 50 

C 25 H 52 

C 26 H 54 

O 

00 

CO 

X 

0 


Mol Wt 

198 

212 

226 

240 

254 

268 

282 

296 

310 

324 

338 

352 

366 

390 

amu 

A|_gH 

72 

77 

82 

86 

91 

96 

102 

107 

112 

117 

122 

127 

132 

110 

kJ mol' 1 

A s-gH 

118 

109 

135 

126 

153 

143 

170 

170 

189 

192 

208 

211 

224 

150 

kJ mol' 1 

A S .|H 

46 

32 

53 

40 

62 

47 

68 

63 

77 

75 

86 

84 

92 

40 

kJ mol' 1 

ln(P v (T r )) 

0.46 

-0.54 

-1.66 

-2.77 

-3.91 

-5.04 

-6.17 

-7.31 

-8.45 

-9.55 

-10.65 

-11.72 

-12.78 

-11.50 

ln[Pa] 

Pv(T r ) l-g 

1.6E+00 

5.8E-01 

1.9E-01 

6.3E-02 

2.0E-02 

6.5E-03 

2.1E-03 

6.7E-04 

2.1E-04 

7. IE-05 

2.4E-05 

8.1E-06 

2.8E-06 

1.0E-05 

Pa 

T r 

298 

298 

298 

298 

298 

298 

298 

298 

298 

298 

298 

298 

298 

293 

K 

T b 

527 

544 

560 

575 

589 

603 

617 

630 

642 

653 

664 

675 

664 

660 

K 

T m 

279 

283 

291 

295 

301 

304 

309 

314 

316 

321 

324 

327 

330 

223 

K 


Vaporization or sublimation at +20 C relative to that at -60 C 


D t p v (293K) 

3.5E+02 

1.3E+02 

4.2E+01 

1.2E+01 

2.2E+00 

7.0E-01 

1.0E-01 

2.5E-02 

4.5E-03 

1.0E-03 

1.6E-04 

4.5E-05 

8.0E-06 

5. IE-03 

pg cm' z s' 1 

D t p v (21 3K) 

1.3E-05 

1.2E-05 

5.2E-08 

5.3E-08 

1.5E-10 

2.2E-10 

5.1E-13 

1.2E-13 

1.2E-15 

1.6E-16 

2.3E-18 

3.9E-19 

9.5E-21 

9.4E-11 

pg cm' 2 s' 1 

D t p v ratio 

2.6E+07 

1.1E+07 

8.0E+08 

2.3E+08 

1.5E+10 

3.2E+09 

2.1E+11 

2.1E+11 

3.9E+12 

6.1E+12 

7.2E+13 

1.1E+14 

8.5E+14 

5.4E+07 


P v (293K) 

9.6E-01 

3.4E-01 

1. IE-01 

3. IE-02 

5.4E-03 

1.7E-03 

2.4E-04 

5.7E-05 

9.9E-06 

2.2E-06 

3.5E-07 

9.4E-08 

1.6E-08 

1.0E-05 

Pa 

P v (21 3K) 

3.1E-08 

2.7E-08 

1. IE-10 

1. IE-10 

3.1E-13 

4.4E-13 

1.0E-15 

2.4E-16 

2.2E-18 

3.0E-19 

4. IE-21 

7.0E-22 

1.7E-23 

1.6E-13 

Pa 

P v ratio 

3.1E+07 

1.2E+07 

9.4E+08 

2.7E+08 

1.8E+10 

3.8E+09 

2.4E+11 

2.4E+11 

4.5E+12 

7.2E+12 

8.5E+13 

1.3E+14 

1.0E+15 

6.4E+07 


ln(P v ) delta 

17.24 

16.34 

20.66 

19.43 

23.59 

22.05 

26.21 

26.21 

29.14 

29.60 

32.07 

32.53 

34.53 

17.97 


A c . g H eff. 

112 

106 

134 

126 

153 

143 

170 

170 

189 

192 

208 

211 

224 

117 

kJ mol' 1 



Simulations assume the 
known temperature 
dependence of one of several 

organic materials; p(To) 
(normalization) is a free 
parameter 
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Optical Bench 
Closeout 


Geometry Model 


Collimator 



Camera Body 

(ACIS Housing) 


From NeilTice/LMC 

Thermal Desktop: 

RADCAD to calculate geometric 
view factors 

SINDA FLUINT to calculate 
temperatures 

738 nodes, 121 OBF-1,203 OBF-S 


OBF-I & S 

i 



ACIS Door 
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Geometry Model 



.000 


REVISIONS 

REV. 

ECO NO. 

DISCRETION 

CHECKED 

APPROVED 

DATE 

B 

36-122 

SEE ECO 

PCI 

RFG 

2/6/95 

C 

36-335 

SEE ECO 

PCT 

RFG 

10/23/95 

D 

36-436 

SEE ECO 





1. MACHINED SURFACES 125/ MAX. EXCEPT 
AS NOTED. 

2. BREAK SHARP EDGES .01 RADIUS. 


4. HECTROLESS NICKEL PLATE PER MIL-C-26074, REV. D, 
CLASS 1 FOLLOWED BY GOLD PLATE PER MIL-G-45Q24 
TYPE III, CLASS 3, GRADE A. 

5. IRIND NO’s 1 AND 4 ARE TO BE PRECISION CLEANED 

TO LEVEL 100A PER MIL- STD-1 246, AT LMA, PRIOR TO 

6/25/2OWLi2A3i^r20)FPMN PDFWriter(D).pc 3 




E?3 


INTERLOCKING DETAIL 


STORAGE CONTAINER 36-60202 IS TO BE PRECISION CLEANED 
TO LEVEL 100A PER MIL-STD-1246 PRIOR TO INSTALLATION OF 
ASSEMBLY 36-10116 FOR SHIPMENT TO MIT. iq 

8£^ 

82 


DETAIL C 

INTERLOCKING DETAIL 
SCALE: 4/1 


NAME 

DATE 

DRAWN J.COOK 

8/29/94 

check edp.c.TAPPAN 

1/23/95 

aprroved^pqoeice 

2/6/95 

RELEASED K _ J| gOEJJg 

2/6/95 

VDEieHT 


80230 | 

a 2/11 1 



IMAGER 

OR- ASSEMBLY 


36-10116 
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Thermal Model 


T~ -62 °C 



T[°C] 

0 

- 7.2 

- 14.4 

- 21.6 

- 28.8 

- 36.0 

- 43.2 

- 50.4 

- 57.6 

- 64.8 

- 72.0 


T~ -69 °C 



AH heaters ON (t<9 yrs) AH heaters OFF (t>9 yrs) 
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Thermal Model 


AH heaters ON AH heaters OFF 



8 = 0.40 



8 = 0.05 


-40 °C 
-50 °C 

-60 °C 

-70 °C 
-80 °C 
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Chandra Contaminant Migration Model 

Results 

Mass Column of Octadecane (C 18 H 38 ) at t=9 years 



low volatility (0. 10) 

“view-factor” dominated: central 
regions have highest accumulation 
because center views more nearby 
cold surfaces, pattern is asymmetric 



high volatility (2.50) 

“thermal” dominated: warm central 
regions begin to clean, pattern follows 
local temperature distribution with 
more material near cold edges 
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Chandra Contaminant Migration Model 

Results 

Mass Column of Octadecane (C18H38) at t =9 years 

As volatility (mass 
vaporization rate) 
increases, center begins 
to clean. At very high 
rates, even accumulation 
at edges is low. 

Observed edge/center 
differences reproduced 
only for a small range of 
volatilities (here, I to 3 
times CisHbs) 

0.1 1 10 100 1000 10 4 

Vaporization Rate (scaled from Dioctyl Phthalate) 
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Chandra Contaminant Migration Model 

Results 



Source rate and material 
properties chosen to 
reproduce AV’s center 
accumulation up to AH 
heaters OFF (t<9yr) 

Model fails to reproduce 
observed edge (and hence A) 
particularly after heaters 
turned off 

Need a second contaminant 

* higher volatility 

* accumulates mainly during 
colder heaters-off times 
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Chandra Contaminant Migration Model 

Results 


A second contaminant can be 
found that accumulates 
mainly on the edges. 

However, time evolution does 
not follow observations well 
and post AH-heaters-OFF 
still problematic 

* low 8 at early times? 

* different source rate time 
dependence for 2 nd material? 


0 2 4 6 8 10 12 14 

Time since 2000 (years) 
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Chandra Contaminant Migration Model 

SUMMARY 

♦ High volatility cleans OBFs and low volatility produces a high 
build-up at OBF centers; only a narrow (factor of 2 or less) 
volatility range produces the observed spatial pattern 

♦ Simulations predict less accumulation above outer S-array CCDs; 
this may explain, in part, gratings/imaging C/MnL discrepancies 

♦ Simulations produce a change in center accumulation due solely 
to DH heater ON/OFF temperature change; but a 2 nd contaminant 
and perhaps a change in source rate is also required 

♦ Emissivity s may depend on thickness; another model parameter 

♦ Additional physics, e.g., surface migration, is not warranted at this 
time 

♦ At t — 14 yrs, model produced 0.22 grams of contaminant, 0.085 
grams remaining within ACTS cavity; 7% (6mg) on OBFs. 
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